Weak measurement: Effect of the detector dynamics 
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The contribution of the detector dynamics to the weak measurement is analysed. According to 
the usual theory [Y. Aharonov, D. Z. Albert, and L. Vaidman, Phys. Rev. Lett. 60, 1351 (1988)] 
the outcome of a weak measurement with preselection and postselection can be expressed as the 
real part of a complex number: the weak value. By accounting for the Hamiltonian evolution of 
the detector, here we find that there is a contribution proportional to the imaginary part of the 
weak value to the outcome of the weak measurement. This is due to the coherence of the probe 
being essential for the concept of complex weak value to be meaningful. As a particular example, 
we consider the measurement of a spin component and find that the contribution of the imaginary 
part of the weak value is sizeable. 



I. INTRODUCTION 

The concept of weak value was introduced in [l[ . It is 
the complex number A w = (Sf\A\Si)/(Sf\Si) in terms of 
which one can express (A), the average result of a mea- 
surement of an observable A preceded by a preparation in 
the state \Si) and followed by a postselection of the state 
\Sf), provided that the interaction between the system 
and the detector, which we shall call the probe as a re- 
minder of its quantum nature, is weak enough compared 
to the coherence scale of the latter [2] . 

Under the assumptions of for a weak analogue of 
an ideal von Neumann measurement, the average value is 
given by (A) = Re(A w ). A surprising result is that this 
average value can lie well outside the range of the eigen- 
values of A. This fact has been confirmed experimentally 
@: EL in optics. Also, the formalism of the weak value 
was proved to describe some relevant phenomena in tele- 
com fibers |6], and to be connected with the response 
function of a system Q. The possibility of performing 
a weak measurement in solid state systems is currently 
being investigated [|[ . In Ref. [l[ the initial state of the 
probe is assumed a pure gaussian state, with a special 
choice of the phase, and the free evolution of the probe 
is neglected. Since the coherence of the probe is an es- 
sential requisite for the weak value to be significant, and 
since the Hamiltonian evolution induces a relative phase 
between different components of the state of the probe, 
the latter assumption seems unrealistic, especially for a 
measurement lasting a finite time. 

In this paper we calculate (A) for any initial state of 
the probe and for any interaction strength [Eqs. 
In the limit of a weak interaction, we show that including 
the free evolution of the probe gives rise to a contribu- 
tion oc Im(A w ) to (A) [Eq. JSJ]; this, generally, does 
not change the main property of the weak measurement, 
namely that (A) can lie well outside the spectrum of A. 
We then consider, as a special example, a probe prepared 
in a general gaussian state, including the state assumed 
in [l[ as a particular case, and we provide additionally an 
expression for the variance (AA 2 ), Eq. (fT5|) . Finally, we 



take A to be a spin component, as a simple illustration; 
we discuss the regime where the weak value does not ap- 
ply, providing formulas for the extrema of (A), (AA 2 ) as 
a function of the postselection [Eqs. (|17I18I20I21[) ] . 



II. MEASUREMENT STATISTICS WITH PRE 
AND POSTSELECTED STATES 

Let us consider a quantum system prepared, at time ti, 
in a pure state \Si) (preselection). We denote the Hamil- 
tonian of the system by H sys . The system interacts, at 
time Ti > ti, with another quantum system, the probe, 
through Hi n t — ~g(t)^qA, where A is an operator on 
the system's Hilbert space, q on the probe's, and g(t) 
is a function vanishing outside a finite interval [Ti,Tf], 

with J g(t)dt = 1. For the measurement to be ideal, if A 
is not conserved, the interaction must be instantaneous, 
Tf = Ti, g(t) — S(t — Ti); otherwise, if A is conserved, 

i.e. A, H sys = 0, the interaction can last a finite time, 

and the measurement is a non-demolition one 9]. The 
probe is prepared, at time t p < Ti, in a state described 
by the density matrix p, and its free evolution is gov- 
erned by H p (p), where p is the conjugate observable of 
q. The operator p is the observable of the probe that 
carries information about the measured quantity A. We 
notice that, in order for the measurement to be ideal, p 
must be conserved during the free evolution of the probe, 
and change only due to the interaction with the observed 
system. At time tf > Tf, a sharp measurement ff2j of an 
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FIG. 1: A schematic view of the measurement with pre- and 
post-selection, the horizontal direction representing increasing 
time. 
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observable Sf of the system is made, giving an outcome 
S, corresponding to the eigenstate IS"). At time t > Tf 
a sharp measurement of the observable p is made on the 
probe. Since p is conserved during the free evolution of 
the probe, this value will not depend on the time t. The 
observed value of non-conserved quantities, by contrast, 
would depend on t [HI . Finally, only those trials in which 
the last measurement on the system gave an arbitrarily 
fixed outcome S — Sf will be selected (postselection). 
The procedure detailed above describes a measurement 
with pre- and post-selection. In Fig.© we provide a 
sketch of the procedure. 

The joint probability of observing the outcome p for 
the probe, at any time t > Tf and Sf for the system, at 
time tf, is given by Born's rule 

p(p,s f \p,s i )= [ dp dp' Q A H ^- H ^ T ^' h 



x p(p ,p' Q )(S f \{p\U Ti ,T f \Si) bo) {SWSAuTj \s f ) \p) 

(1) 

where U is the time evolution operator for 
H sys + Hp + Hi n t, we introduced p{p,p') the 
probe density matrix in the \p) basis, and 

\Sij) := exp $y-iH sys (Tij - Uj)/h^\Sij). After 

introducing twice the identity resolved in terms of the 
eigenstates of A, we obtain 

T( P , S f \p, S^ = J>(p - Xa >P~ Aa')e-^"-*«'] 

a, a/ 

x (SflaXalSMSila'Ha'lSf), (2) 
with a, | a) the eigenvalues and eigenvectors of A and 

$a := \ I ^ H p(p- Xa J s f ds '9{s')^j ■ (3) 

For an instantaneous interaction, Tf in Eq. ([3]) and else- 
where in this paper should be interpreted as being a time 
infinitcsimally later than the interaction. In deriving 
Eq. ©, we exploited 

(p|re-^o d ^^ s ^^|po) = s(p-P0- ^ dsf(s) 

x e - i Io dsH P^Sl ds ' f( s ' ))/ h ^ (4) 

We notice that if no postselection were made (i.e. if one 
would sum over the final states \Sf)), the off-diagonal 
elements of p(p,p') would not contribute to Eq. ©. 

The conditional probability of obtaining outcome p, 
given that the state has been postselected in Sf, is 



P(p\S f ,p,Si) = 



r(p,s f \p,Si) 
fdp'vd/^fip^iY 



(5) 



where we applied Bayes' rule, and the expected value 
inferred for the system through observation of the probe 

[ } - A jdp'V( P ',Sf\p,S t ) ■ [ > 



Since what matters is the deviation of the pointer p 
from its unperturbed expected value J dp pp{p,p), we 
can set the latter to be zero without loss of general- 
ity. We notice that the inference assigning the quantity 
A = p/X to the system when observing the probe to have 
the value p is valid only when initially the probe has a 
precise enough value of p close to zero. We can, how- 
ever, keep assigning the value to the system even when 
the probe is not sharply prepared around p = 0, and 
say that we observed a value A for the measured sys- 
tem (and correspondingly we shall define a probability 
V(A) := XP{p/\)). This value is not necessarily one of 
the eigenvalues of A, and, as shown in it can even 
lie outside the range [a m i n , o-max] (for this reason we are 
indicating with a the eigenvalues of A and with A the 
outcome of each individual measurement). 



III. WEAK MEASUREMENT 

A measurement is weak when the coupling A is small 
compared to the coherence length of the probe, i.e. to 
the range of \p — p'\ within which p{p,p') vanishes. This 
can be evinced from Eq. ([2]). 

In the following, we shall assume that p(p,p') is ana- 
lytic in a neighborhood of p = p' . Then, to lowest order 
in A, the denominator in Eq.fJS]) is J dpV(p, S f\p, Si) = 
\{Sf\Si)\^ . Before analysing the numerator in Eq.©, we 
rewrite p(p,p') — F(p,p') exp {ia(p,p')}, with F, a sym- 
metric and antisymmetric real functions, respectively. 
We have then that the numerator in Eq.© is 

dpfp(p,5 / |p ! S i )~-^<5 / |a)<a|5 i )(5 i |a')<a'|5 / ) 



=Re 



3 ((5 / |i|^)<^|5 / ))~pG(p)/m((^ / |i|5 J )<5 J :|5 / )), 

where V(p) := F(p,p) is the initial distribution of the p 
observable of the probe, 



G( P ) 



2At dH p (p) 2 da(p,p>) 



dp 



dp 



p'=p 



ds l ds'g(s') , 



(7) 



and the bar symbol denotes the average over V(p). After 
introducing the weak value, i.e. the complex number 
A w := {S f\A\Si) / (S f\Si) , the average value is 



(A) ~ (A) := Re{A w } - pG(p) Im{A w } . (8) 

Eq.© holds as far as the product of the prepared and 
the postselected state is larger than the first nonvanishing 
contribution in the A expansion for the denominator. In 
the latter case, one should keep the latter contribution 
as well. 
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IV. COMPARISON WITH PREVIOUS RESULTS 

We notice that the contribution of the imaginary part 
has been generally overlooked in the literature, due to 
the neglecting of the Hamiltonian of the probe and to 
the choice of a very special phase a(p,p') = 0. On the 
other hand, it has been proved [H,[I3| that observing the 
q variable of the probe one gets an average value which 
is proportional to the imaginary part of A w . This is true 
only if one neglects the time evolution of the probe from 
preparation to observation. When this evolution is ac- 
counted for, the observed value of q depends also on the 
details of the free Hamiltonian of the probe and on the 
time of observation. To the best of our knowledge, the 
first paper to point out that Im(A w ) contributes to (A) 
was reference [111 ]. There, however, the readout variable 
p (which in the notation of [llj is actually q) is not con- 
served during the free evolution of the probe. Thus the 
results presented in [ll[ hold only if the system-probe 
interaction is instantaneous, and if the probe is read im- 
mediately after the interaction. Indeed, if the interaction 
is instantaneous, the probe is prepared at time t p and it is 
observed at time t, the central result, Eq.(17), of Ref. [ll[ 
should be substituted by (we use our notation p <-» qj 



(M) ~ (M(t)) p + iX 



ReA w ( M(t),q(Ti) ) p 



ImA, 



({{M{t)A{Ti)}) P -2{M{t))MTi))p 



(9) 



where [•,•] ({•,•}) denotes (anti)commutator, (M) p = 
Tr{pM} is the trace taken with the probe density matrix 
at time t p , and M(t) is the probe observable M evolved 
with H p in the interval [£ p ,i]. Generally, with a probe 
Hamiltonian H p = p 2 /2M p + V(q), one can no longer link 
the contribution of Im{A w } to (p) with the derivative of 
the variance of q, unless t = T^. However, if the probe 



V. PROBE PREPARED IN A MIXED 
GAUSSIAN STATE 

So far, in the literature on the weak measurement, the 
probe was assumed to be prepared in a pure gaussian 
state. The corresponding density matrix p{p,p') is char- 
acterized by the identity between the scale in \p — p'\ 
over which its off-diagonal elements vanish (the coher- 
ence length scale) and the scale over which the diagonal 
elements decay going away from the zero value (the clas- 
sical uncertainty spread in p). We shall consider a more 
general gaussian distribution 

p-{(p+p') 2 /8AP 2 +(p-p') 2 /8<5p 2 -i(p-p')/ 2 P4 

P(P,P') = - = ■ (12) 



2ttAP 

Here AP is the initial spread and 8p the coherence scale. 
Positive semidefineteness requires Sp < AP. We assumed 
a phase linear in p, with p$ a scale. The linear phase 
such chosen defines the center of the Wigner function in 
the coordinate Q, Qq :— h/2p l p. We take a quadratic 
Hamiltonian for the free probe H p (p) = p 2 /2M p , and 
we define a further scale pn '■= \J hM p /2At, with At 
defined by Eq.([7]). We stress that the presence of Vh 
can make this scale the smallest one. We have then that 
the average det ected v alue of the observable A is given 
by Eo, © with pG(p) = AP 2 /p 2 H := k 2 . The results of 
Ref. m are recovered for AP = Sp, and ph,P4> — ► oo. We 
also provide an expression for the variance of A 



(AA 2 ) 



AP 2 



1 



1 



ReAA 2 



K 2 ImAA 2 w , 



(13) 



where we introduced AA 2 ^ :— (A 2 ) w — (A w ) 2 , with 
(A 2 ) w := (Sf\A 2 \Si)/{Sf\Si). We notice that there is 
always a large contribution AP 2 /X 2 , due to the ini- 
tial spread in p. The calculated variance differs from 
Eqs. (24,25) of Ref. p3: even for k = (which is 



Hamiltonian is H p = p 2 /2M p , we have that p(t) = p{Ti) the limit considered in [lfj), Eq. (filj|) allows (AA 2 ) < 

AP 2 /\ 2 . We also stress that if the last two terms in 
Eq. lfTB")) take a small value, this does not imply that a 
single measurement can reveal the weak value: A is in- 
ferred from the observed p of the probe through A = p/X; 
since p has a spread of order AP, the value of A observed 
in each individual measurement will vary with a spread 
of order AP/A, which is large by hypothesis. 



and thus 

M 

(p) ~ XRe{A w } - X-ZpIm{A w }, (10) 
n 

where (3 := dVarq(t)/dt\ t=T . The formula agrees with 
the more general Eq. ©, since for a generic H p (p) 



dVarq(t) 
Jt 



H^M G (p) 



dp 



dp 



which reduces to dVarq(t)/dt — hpG(p)/M p . For a 
generic Hamiltonian H p (p), if instead of observing p, 
one observes the "velocity" operator V — V(p) :— 
dH p (p) I dp, one has 



(V)~(V) p + X& p Re{A w }~X(3Im{A w }. 
dp 



(11) 



VI. ILLUSTRATION: WEAK MEASUREMENT 
OF A SPIN COMPONENT 

As a specific example, we consider a measurement 
of spin components. We assume that the interaction 
between the spin and the probe lasts a finite time T, 
that g(t) = 1/T during the interaction and zero other- 
wise, and that the probe is prepared in the state p of 
Eq. (|12p immediately before the beginning of the inter- 
action. Then Eq. © gives At = T/2. We take the spin 
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FIG. 2: (A) as a function of 7 for 8 = ir/2, and <j> = 7r/2 (top 
figure), (j> — <f>* — 7r/4 (bottom figure). (The inset shows how 
the angles were defined.) The full line depicts the exact value, 
the dotted one depicts the approximate value including the 
contribution of Im{A w }, and the dashed line the contribution 
of Re{A w }. The thin vertical and horizontal lines correspond 
to 7* = 7T - X/u, (A) m = e(4>) v/X [see Eqs. (|17I18|| ] . In 
the range 7 = [0,3], the dotted and the full line practically 
coincide. We assumed X/Sp — 0.01, X/p^ — 0, 8p — AP — pn- 



to have been preselected in the state up along direction 
rij and postselected in the state up along a direction n/, 
while A = n ■ er. Then we obtain from Eqs.foT 



tr=±l I 



V2ttAPN 
(711-110(1 + an . n/ ) e -(f- ffA ) 2 / 2Ap2 



(nx rij) • (nxn/) cos [AG(p)] 



n • (iijXn/) sin [AG(p)] 



; -p 2 /2AP 2 -A 2 /25p 2 I ^ 14 ) 



N = 1 + n 4 • n/ + e _A /2l/ sin ( — )n • (n» x n/) 
A 



A 



1 - e 



(n x n 4 ) • (n x n/), (15) 



with 1/ := [1/V + k 4 /AP 2 ]~ 1/2 < AP/yTT^. 
The exact average value is 



(A) =< n-n.i + n-Uf — k 



_ K 2-\ 2 /2,y 2 



cos 



sin I — I (nxn,*) ■ (nxn/) 



(-]n-( n ' xn /) 

w 

/ N. (16) 



To lowest order in A, (A) is given by Eq. §8§ with A w = 
n- [n, + n.f + ixii x n/] / [1 + n^ • n/]. Interestingly, when 
n lies in the plane orthogonal to the bisector of n, and 
iif, A w is purely imaginary. This setting of the weak 
measurement can hence be a testing ground to detect 
the contribution of the imaginary part. 




FIG. 3: (a) The variance of A as a function of 7 for fixed 
6 = 7r/2, <f) — 7r /4- The dashed vertical and horizontal lines 
correspond to the positions and the values of the extrema 
[Eqs. (120121)) ]. (b) The probability distribution for 6 = ir/2, 
<f> = 7r/4 for some significant values of 7. The parameters are 
the same as those of Fig. @. 



Without loss of generality, we take n^ , n to dehne the 
XZ plane, with the former as the Z-axis, and the latter 
forming an angle 9 € [0, ir}. The direction n/ is define d by 
the azimuthal and polar angles 7 € [0, 2tt], <j> € [0, 7r] [Tj]. 
Then A w = cos 9 + sinfle - ^ tan (7/2). In Figure [U we 
plotted (A) as a function of 7 for fixed 9, <j>. We compare 
the exact value [Eqs. (jl6H5|) ] . the approximate value (A) 
given by Eq. ©, and Re(A w ). For fixed 9 ^ 0, n, and <j>, 
(A) reaches extremal values for 7* =71 — 77*, 



smf 




(17) 



with e(4>) := cos <fi + k 2 sin <f>. The extremal value is 



(A) m = 



cos(/>£'(0) ± \ e(4>) 2 , (l + /t 4 )(cos0) z 



(18) 



with the prime meaning differentiation. Generally, the 
upper sign solution in Eqs. (I17I18[) holds as far as e((f>) ^> 
A/p^,, X/v, in which case the extremum is no longer found 
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close to 7r, but 



sin 7 



-2K 2 A s in^ E (^)+ e '(0)Acose) 



(k 2 ± sin 6» 



e($+s'(0)£cc*0 



(19) 

and (A) m ~ cos 6*, while the lower sign solution converges 
to 77* ~ —XsmO/p^y/l + k 4 , 



<^>t 



2K 2 X/p <j> 



k 4 A 2 /(1 + k 4 )p 2 +A 2 /V 2 



There are two exceptions to this: (i) For ^> v, 
if ~ ±wa.0X/v, and (A) OT = ±e(4>)v/X. (ii) For 

k 2 = 2AP 2 At/hM p < 1, 77* ~ ±sm.0XJl/v 2 + 1/p 2 



— sin 



and (A) m =cos(0)/A [±^1/^ + i/ p 2 

The extremal value for (A) as a function of both 7, <f> 
has an involved expression, except for p^, > y, when 
the location of the extremum is 7* = n =p &m.0X/v, </>* = 
Arctan(K 2 ), and (A) m = ±\/l + n A i>/\. In the same 
limit, we have that the minimum of the spread is reached 
for 77* ~ ±\/3sm9X/v, 4>* = Arctan(n 2 ) 



(AA 2 



AP 2 -(1 + K 4 )^ 2 /4 > 3AP 2 ^ 



Its maximum is reached for 7* = 7T, and it is 
AP 2 + 2(1 + k> 2 AP 2 



(AA 2 



A 2 



< 3- 



A 2 



(21) 



A 2 



4 A 2 



We plot the probability distribution for three values of 
7: close to 7 = 7r the distribution has two peaks, each of 
order 100 for the choice of parameters made. While the 
average value goes to zero when 7 gets very close to 7r (n— 
7 <C X/v), the probability density of observing a value 
in the range [—1,1] is rather small: in each individual 
measurement, it is likely that the value of \A\ will be 
much larger than unity. 



VII. CONCLUSIONS 

We have showed that accounting for the dynamics of 
the probe in the weak measurement leads to an observ- 
able deviation of the average value from the real part of 
the complex weak value defined in Ref . |l| . We have also 
derived an expression for the spread, and, in the case of 
spin, we have individuated the locations and values of 
the extrema of (A) . 
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